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Abstract 



In this paper we summarize the minimal supersymmetric standard model 
as well as the renormalization group equations of its parameters. We proceed 
to examine the feasability of the model when the breaking of supersymmetry 
is parametrized by the soft terms suggested by supergravity theories. In 
such models, the electroweak symmetry is exact at tree level and is broken 
spontaneously at one loop order. We make the additional assumption that the 
GUT-inspired relation mf, = ttIt- be valid at the scale where the gauge coupling 
constants unify, which constrains the value of the top quark mass. For all types 
of soft breaking terms expected in supergravity theories, we present the results 
of numerical runs which yield electroweak breaking at the required scale. 
These yield not only the allowed ranges for the soft supersymmetry breaking 
parameters, but also the value of the supersymmetric partner' masses. For 
example in the strict no-scale model, in which global supersymmetry breaking 
arises solely from soft gaugino masses, we find that Mt can be no heavier than 
~ 127 GeV. 
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I. INTRODUCTION 



In the last few years, it has become apparent, using the ever increasing accuracy in the 
measurement of the strong couphng, that Supersymmetry (SUSY) affords an elegant means 
to achieve gauge coupling unification at scales consistent with Grand Unified Theories 
(GUTs) Whereas in the Standard Model (SM) the three gauge couplings unify "two 

by two" forming the "GUT triangle," in the simplest Minimal SUSY Extension of the SM 
(MSSM) these gauge couplings spectacularly unify at a point (within the experimental errors 
in their values). Given that the scale of unification in these models is generally above the 
lower bound set by proton decay, these supersymmetric grand unified theories (SUSY-GUTs) 
have regained increasing interest. Constraints from the Yukawa sectors of such models have 
also yielded interesting predictions for various low energy parameters including the top quark 
mass IHH]. 

The analyses mentioned above have employed the Renormalization Group (RG). In the 



first part of this two part series [|T3|, we reviewed the use of the RG for the purposes stated 
above in the framework of the SM. We discussed the aspects of data extraction involved 
in arriving at boundary conditions for the renormalization group equations (RGEs). In 
attempting to make the work as complete as possible, we included all RG /^-functions to 
2-loops making no approximations in the Yukawa sector. Plots of the evolution of all the 
SM parameters to this order were presented. Furthermore, we discussed threshold effects in 
the running fermion masses, an analysis which to our knowledge had not been adequately 
treated in the literature. 

In the present paper we generalize this analysis to the minimal supersymmetric extension 
of the standard model. We include the two loop renormalization group equations for the 
parameters of the model. We only consider the case where supersymmetry is broken by 
soft terms in the Lagrangean, and use the renormalization group equations for these soft 
breaking parameters at the one loop level. One of the purposes of this paper is to determine 
the range of these parameters which leads to electroweak breaking with the correct value of 
the Z-mass. We discuss in some detail the aspects of the effective one loop potential which 
are relevant to this breaking. The mass formulas for the sfermions, Higgses, charginos, 
and neutralinos are also presented for completeness. Various boundary conditions at the 
unification scale in these minimal low energy supergravity models are discussed. We then 
describe our numerical procedure. The treatment of thresholds and the "special" form of 
the /5-functions needed is then discussed . Although similar analyses have appeared in the 



literature |]T^|-0, we feel it important to present our results in a way that directly compares 
all allowed sets of values for the soft supersymmetry breaking parameters. Our purpose is 
to use the results of our analysis to find to which extent the low energy data constrains the 
types of supersymmetry breakings. Supersymmetry presents us with a way to extrapolate 
low energy data to the deep ultraviolet, where they may be compared with the underlying 
theory, be it GUTs and/or Superstrings. 
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II. MINIMAL SUPERSYMMETRIC STANDARD MODEL 



In the minimal supersymmetric extension of the standard model, every particle has a 
supersymmetric partner, their spins differing by a half |jl8[. In addition, two Higgs fields 



with opposite hypercharges are needed. Since the superpotential cannot consist of fields 
and their complex conjugates, phenomenology requires two Higgs fields to give mass to the 
charge +2/3 and —1/3, —1 sectors, respectively. Remarkably, two Higgs fields are also 
needed for both chiral and SU{2) global anomaly cancellation. 

For renormalizable theories, the superpotential can have at most degree three interac- 
tions. The superpotential for the MSSM is (suppressing the SU (2) and Weyl metrics) 

W = flYjuQ + lYd^dQ + tY,<l,L , (2.1) 

where the hat indicates a chiral superfield and the overline denotes a left handed CP con- 
jugate of a right handed field, ip = ia2'ip*ji. The usual Yukawa interactions are accompanied 
by new Yukawa interactions among the scalar quarks and leptons and the Higgsinos in the 
supersymmetric Lagrangean. There are also new gauge Yukawa interactions involving the 
gauginos. The new purely scalar interactions form the scalar potential which is positive 
definite in supersymmetric theories. The scalar potential will be discussed in a subsequent 
section. A remarkable aspect of supersymmetry is that all these new interactions require no 
new couplings. Table I displays the 5't/(3) x SU{2) x f/(l) quantum numbers of the chiral 
(all left handed) and vector superfields of the MSSM. 

This superpotential exhibits two anomalous chiral global symmetries. The first, called R- 
symmetry, is characteristic of supersymmetric theories with cubic superpotentials; it implies 
massless gauginos, it must be broken; this is simply done by adding gaugino masses, breaking 
it down to a discrete R-parity. The other is an exact Peccei-Quinn (PQ) symmetry which, 
with this superpotential, can only be spontaneously broken at the electroweak scale. This 
is well known to lead to an axion, heavy enough to be ruled out by experiment. Hence the 
superpotential has to be improved so as to avoid this difficulty. This can be done in several 
ways. 

• The standard way is to break the PQ symmetry explicitly through the addition of the 
following supersymmetric-invariant term to the superpotential 

^JiK^d ■ (2.2) 

This represents the simplest way to avoid the phenomenological disaster just described. It 
is intriguing that, in order to achieve the required electroweak breaking scale, the value of 
yU turns out to be (as we shall see in the next section)of the same order as the soft breaking 
parameters. At this level, it is not possible to explain this coincidence, since these terms 
break different symmetries. 

• Alternatively, the /i can be interpreted dynamically as the vacuum expectation value 
of a singlet chiral superfield, iV, through the following interactions which also preserve 
supersymmetry and break PQ explicitly in the Lagrangean 

XN^u^d - XqNNN . (2.3) 
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TABLE I. Quantum Numbers 
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The effective can be identified with A < >. This approach provides a natural explana- 



tion for ~ 0{Mw), if A, Aq ~ 0(^2) @ 



• A third possibihty implements non-electroweak spontaneous PQ symmetry breaking and 
leads to an invisible axion, but it needs a special scale of the order ~ 10^^ GeV. This 
introduces another hierarchy of scales. This model is parametrized by the following addition 
to the superpotential 

AX<I,<I>, + AoXff . (2.4) 

• Finally, one can also break both the PQ symmetry and supersymmetry by adding to the 
Lagrangean the explicit soft breaking term, 

m^$„$rf . (2.5) 

This term can be put in by hand or generated in the low energy supergravity model for 



sufficiently general couplings of supergravity to the standard model ||20|| . 

If the superpotential contains the /i term, then this soft term will be generated even for 
the case of minimal coupling. In such a case, it is convenient to set m| = S/x, where B is 
the soft bilinear coefficient. In the following we will adopt the simpler scheme of the /i term, 
leaving the other possibilities to future investigations. 



III. MINIMAL LOW ENERGY SUPERGRAVITY MODEL 

Since no super particles have been observed experimentally, supersymmetry, if truly 
present in nature, must be broken. One way to accomplish this breaking is to add to the 
Lagrangean soft breaking terms. 

The general soft symmetry breaking potential for the MSSM can be written (including 
gaugino mass terms) 

Vsoft = ml^^i^^ + m|^<l>;^<l>^ + 5/i($„$d + h.c.) 

i * 

+ E( ^iyu'^^^uQ, + A'M%^aQj + + h.c. ) , (3.1a) 
1 3 

Vgaugino = T^Yl ^l^l^l + ^"C- , (3.1b) 
^ 1=1 
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where Vgaugino is the Majorana mass terms for the gaugino fields, A; (suppressing the group 
index), corresponding to U{1), SU{2), and SU{3), respectively. 

We see that there can be sixty-three different soft symmetry breaking parameters. In 
addition, these could each be introduced in the effective Lagrangean at their own scale! 
Thus the supersymmetry breaking section runs the risk of being more complicated than the 
standard model itself, which only has a mere eighteen parameters! 

It is clear that we need a further organizing principle to describe supersymmetry breaking. 
One attractive possibility is suggested by the coupling of the = 1 standard model to 

= 1 supergravity (SUGRA). The idea is to have supersymmetry broken in a sector that is 
uncoupled to the fields of the standard model, except through the universal interactions with 
supergravity. Then the net result of this picture is to produce in the effective Lagrangean 
of the minimal low energy supergravity model, a specific pattern of induced soft breaking 
terms of the above form, but with far fewer parameters. 

Let us summarize the basic facts of supersymmetry breaking. From the supersymmetry 
algebra, one deduces that spontaneous symmetry breaking occurs if and only if the vacuum 
energy is not zero. In global supersymmetric theories, the scalar potential is a sum of F- 
and Z)-terms. Supersymmetry is spontaneously broken if either the vacuum values of the 
F-term |PT| or D-term are non-zero. A consequence of the spontaneous supersymmetry 



breaking is a massless fermion in analogy with the breaking of an ordinary global symmetry. 
Spontaneous supersymmetry breaking in supergravity occurs via the super-Higgs mecha- 
nism. The Nambu-Goldstone fermion, or Goldstino, associated with the breaking of global 
supersymmetry is eaten by the gravitino thereby providing it with a mass. We will assume 
that the spontaneous breaking of the local N=l supersymmetry is communicated to the 
"visible" sector by weak gravitational interactions from some "hidden" sector. This type of 
spontaneous breaking of supergravity manifests itself at low energy as explicit soft breaking 
terms of supersymmetry. 

The SUGRA Lagrangean is characterized by two arbitrary functions of the fields, a real 
function K (the Kahler potential) that determines the kinetic terms of the chiral superfields, 
and an analytic function /q^, transforming as the symmetric product of the adjoint repre- 
sentation of the gauge group, that determines the kinetic terms of the gauge fields. In terms 
of these, the supergravity scalar potential is 



VsUGRA — e 



K/M 



'dW 1 dK^^ \ fdW* 1 dK^^\ d^K \W\^ 



M2 dchi / \ dd)* M2 d(h* / 86^36^* M"^ 



\f-^pD"D^ , (3.2) 



where M = Mpianck/ V^n and are auxiliary fields. In models with minimal kinetic terms 
for the chiral superfields (fiat Kahler potential, d'^K/d(j)'^d(f)^* = this leads to a common 
(gravitino) mass, mo, for all the scalars of the model. The presence of non-minimal gauge 
kinetic terms {fa/s 7^ Sa/s) implies non-zero masses, M/, for the gauginos at the GUT scale, 
Mx- By further assuming gauge coupling unification, we can take the three gaugino masses 
to be equal. Furthermore, the trilinear soft couplings A^J, A'J, and A^^ are all equal to 
a common value Aq. With minimal chiral kinetic terms, the bilinear soft coupling Bq is 
related to Aq as Bq = Aq — tjiq. This scenario has obvious, desirable features. First, it is 
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very predictive since it has a few parameters accounting for thirty-one new masses. Second, 
the universal nature of the squark and slepton masses at Mx helps to avoid the appearance 
of unwanted flavor changing neutral current (FCNC) effects. In fact, one could argue that 
the absence of FCNCs hints at a universal mass for the scalars. 

In this paper, we will not concern ourselves with the exact nature of the soft breaking 
terms. We will leave to others the issue of finding the exact, more fundamental supergravity 
models which engenders particular sets of vahics. Our philosophy will be to assume the 
presence of these terms and to explore their phenomenological consequences. We shall 
study the effects of varying the values the soft breaking terms over some interesting ranges 
and for some exceptional cases. The scales at which these soft breaking parameters enter 
the effective low energy Lagrangean are determined by theoretical prejudices. In view of 
the unification of gauge couplings with supersymmetry, it seems natural to input the three 
gaugino masses at the GUT scale. We note that the scales at which the sparticle, Higgs 
masses and trilinear couplings enter our Lagrangean could in principle be anywhere between 
the Planck and GUT scales. However, since we are interested in the rough features of the 
models, we have chosen, in the name of simplicity, to enter all the soft breaking parameters 
at the same GUT scale. 



IV. RADIATIVE ELECTROWEAK BREAKING 

The complete scalar tree level potential now appears as 

V ^Vp + VD + Vsoft , (4.1) 
where Vp contains the potential contributions from the F-terms 

+ |^Y„$„ + 5Y^$<i|' + |iY,$rf|2 , (4.2) 
and Vd contains the potential contributions from the D-terms 

= y( ^ QlQ^ - 1^1% + ^f^, - H\U + i% + \^i^u - \^\^d ? 

+|( QlrQi + 1\tU + ^ir^u + ^\T^d f 

+ f ( Ql^Qi - ^lA*^i - t>^% f , (4.3) 

where r = (ti, T2, Ts) are the SU (2) Pauli matrices and A = (Ai, . . . , Ag) are the Gell-Mann 
SU (3) matrices. In general, one must impose constraints on the parameters to avoid charge 
and color breaking minima in the scalar potential. Some necessary constraints have been 
formulated, such as 
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Al < 3(m| + m| + m|J , (4.4a) 
< 3(m| + m? +m|j , (4.4b) 
Al < 3(m| + m| + m| ) . (4.4c) 



However, these relations are in general neither sufficient nor indeed always necessary ^3 
Their derivation involves very specific assumptions about the spontaneous symmetry break- 
ing. 

An appealing feature of the models we are considering is that they can lead to the 
breaking of the electroweak symmetry radiatively [P^-p7|. The one loop effective Higgs 
potential in these models can be expressed as the sum of the tree level potential plus a 
correction coming from the sum of all one loop diagrams with external lines having zero 
momenta 

Vl-loop = VtreeiA) + ^iA) • (4.5) 

The right hand side is A-independent up to one loop. The one loop correction is given by 

^^^(^)=64^^'^^-^'^^"^-2)^ 

1 ..o„ ml 3, 



64,.E(-l)'M2.p+l)m^(ln^--), (4.6) 



p 



where A/1 is the field dependent squared mass matrix of the model and rup is the eigenvalue 
mass of the p*^ particle of spin Sp. The tree level part of the potential is 

VtreeiA) = (A) <|.i( A) ( A) + m2(A)$t (A)<|.,(A) + m^(A)($„(A)$,(A) + h.c.) 
+ ^(*l(A)*„(A)-$i(A)$,(A))^ 



+ ^(<^>I(A)r«f „(A) + $ kA)r$,(A))2 (4.7) 



where 



m2(A) =m|^(A)+/i2(A) , (4.8) 
m^(A)=m|JA)+/x2(A) , (4.9) 

m^(A)=i?(A)MA). (4.10) 

The parameters of the potential are taken as running ones, that is, they vary with scale 
according to the renormalization group. The logarithmic term in the one loop correction is 
necessary in making Vi_;oop(A) independent of A to this order (up to non- field dependent 
terms) . 

Given the low energy scale of electroweak breaking, we must use the renormalization 
group to evolve the parameters of the potential to a convenient scale such as Mz (where the 
experimental values of the gauge couplings are usually cited) thereby making this leading 
log approximation valid. The exact scale is not critical as long as it is in the electroweak 
range. If we define. 
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2 , ^AKi 



^' = ^' + ^' (4-11) 



with Vi = Vd, V2 = Vu and 



then minimization of the potential yields the following two conditions among its parameters 

1 ^ ml-mli^n-p 

2 ^ tan2/5-l ' ^ ' 
where m| = (5?'^ + gl)v'^/2, v'^ = vl + vj, and 

Bfi = -(ml + ml) sin 2P , (4.14) 

where tan/3 = Vu/vd- 

Although results based on the tree level potential cannot always be trusted, one can 
still use it to get some idea under what conditions electroweak breaking occurs. The renor- 
malization group evolution of (see Appendix A) can be such that it turns negative at 
low energies, if the top Yukawa coupling is large enough, whereas m|, runs positive. From 



Eq. (O), the scale at which breaking occurs A^ is set by the condition 

m?(Afc) ml{A,) - m^(Afc) = . (4.15) 

If the free parameters are adjusted properly, then the correct value of the mass [Mz = 
91.17 GeV) can be achieved. 

In the tree level analysis, there is another critical scale A^ that must be considered. It is 
evident from Eq. ([4.7|) that the potential becomes unbounded from below along the equal 
field (neutral component) direction, if 

ml{A,) + ml{A,) < 2ml{As) . (4.16) 



Since m\m\ — m\ > implies ml + m^ > 2m\, condition (|4.16| ) can only occur at scales 



lower than condition ( |4.15| ), so A^ < A?,. From this analysis, one concludes that the tree level 
vacuum expectation values (VEVs) of the scalar fields obtained by minimizing the potential 
are zero above A^, and grow to infinity as one approaches Ag where the potential becomes 
unbounded from below. It follows that the appropriate scale at which to minimize the tree 
level potential and evaluate the VEVs is critical. This scale must be such that the one loop 
corrections of the effective potential may safely be neglected. Only at such a scale can the 
tree level results be trusted. However, there is more than one scale involved, and therefore, 
it is difficult if not impossible to find a scale at which all logarithms may be neglected. 
Indeed, the use of tree level minimization conditions to compute the VEVs at an arbitrary 
scale (e.g., A = Mz) leads to incorrect conclusions about the regions of parameter space 



that yield consistent electroweak breaking scenarios |28|. When AVi is included, however. 



the value of A is not critical as long as it is in the neighborhood of Mz- 
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Reference gives a prescription for arriving at a scale (A) at which the tree level 



and the one loop effective potential results for the VEVs agree. Three qualitatively different 
cases are considered. In Ref . pSf , Msusy parametrizes the superparticle thresholds, then the 



cases can be characterized by the orderings: (a) Msusy < < A;,, (b) A^ < Msusy < ^b, 
and (c) As < Af, < Msusy- In each case, the prescription is to take A = ma.x{MsusY, -^s}- 
Two cases deserve special mention. Case (a) cannot be handled using the tree level analysis 
because Vu,Vd — > oo near A^. Fortunately, phenomenological bounds rule this case out 
anyway. In case (c), there is actually no electroweak breaking. For scales below Msusy, the 
superparticles have decoupled and the effective theory is not supersymmetric. Therefore, 
the running mass parameters of the potential freeze into their values at A = Msusy at 
which scale there is no electroweak breaking. Finally, it must be emphasized that the 
apparent violent behavior of the VEVs with scale in the tree level analysis is an artifact of 
the approximation. The only physical potential is the full effective potential, and it either 
breaks electroweak symmetry or not. If it does, then the scalar fields have non-zero VEVs, 
and these VEVs are non-zero over all scales varying according to the anomalous dimension 
of their respective scalar fields (in the Landau gauge). 

In contrast to the tree level potential, the one loop effective potential is constant against 
the renormalization group to this order around the electroweak scale. The exact scale at 
which to minimize is no longer critical. Moreover the assessment of the masses of the Higgs 
bosons based on the one loop effective potential is more accurate. The tree level restriction 
Mfi < Mz is known not to be valid when one loop corrections, which are large because M^ 
is large, are included in the determination of M^ PP|-PT|. 

In this paper, we do not rely on the tree level analysis, rather we incorporate the one loop 
corrections. We include the dominant contributions from the third family, that is, those of 
the top and stop, bottom and sbottom, and tau and stau p2|j33[]. We choose the Z° mass 



as the scale at which to evaluate the minimization conditions. Equations ( 4.13 ) and ( 4.14 ) 
can be written 



— tan^ /? 1 
tan^ i3 - 1 2 



fi\Mz) = " - -mi , (4.17a) 



2^{Mz) 

where ^ = ^ + (9AVi/(9f^ ^ and used to solve for fi{Mz) and B{Mz) given the value 
of all the relevant parameters at Mz- We note that the form of Eq. ( [4.17a| ) does not fix the 
sign of yU, and a choice for its sign must be made {fi is multiplicatively renormalized; see 
Appendix A). The right hand sides of these equations implicitly involve the VEV at Mz- 
In a consistent scenario it would have the value of v{Mz) = 174.1 GeV. If the parameters 
are such that the right hand side of ( [4.17a|) is negative, then the scenario is inconsistent and 



the electroweak symmetry fails to be broken. 

The validity of using (^4.5|) at Mz hinges on the assumption that there is spontaneous 
symmetry breaking and that Mz = 91.17 GeV. Given a set of values for the input parameters 
(soft terms, etc.) at Mx, we proceed by assuming valid electroweak breaking {i.e., Mz = 
91.17 GeV). The one loop effective potential of Eq. (|4.5|) should then be perturbatively 
valid at Mz with all its running parameters evaluated at this scale. In this way, we are 
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renormalization-group-improving the potential. The vahdity of this assumption is tested by 
the consistency of Eqs. ( |4.17a| ) and ( 4.17b| ). Failure to attain consistency invalidates the 



initial assumption that the given input values at Mx can accommodate a low energy world 
as we know it. 



V. SPARTICLE MASSES 

In the following, we list the tree level mass formulas for the different superpartners. 

A. Sfermion masses 

The mass matrices for scalar matter are constructed from Eq. (|4.1| ). For example, in the 
up squark sector the relevant mass matrix appears as 

\ RiLj RiRj / 

where i,j = 1, 2, 3 are flavor indices and 

MU, = mlS.j + vliYlYj,, - ^{vl - vlWMg" - n^gDS,, , (5.2a) 

Mln, = + vliYlYj,, -^ivj- vl){Y{ua)g'^)5., , (5.2b) 

Ml,L, = ^vaY^^ + (5.2c) 
Mj p =Ml*j . (5.2d) 

Note in Table I that F = Q — T3 in our notation. Similar matrices follow for the other 
sfermions. These mass formulas as well as the ones to follow are given in terms of running 
parameters. The domain of validity of these formulas is at low energies (~ Mz) with the 
parameters taking on their renormalization group evolved values at this scale. 

B. Higgs masses 

If we employ the notation 



then the physical masses of the Higgs at tree level are calculated from the following three 
matrices 



2 9(3?0O)a(5R0O) 2 ^ ^ \^ -1 cot/5; ' 2 ^ \^ -1 tan/5 

^^^3^ = -M,^sm2/5( ^ ^^^A (5.4c) 
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where = m\ + m|, Mjj^ = M\ + rn^r, and = g\v^ j^. The eigenvalues for the first 
matrix are corresponding to the Goldstone boson and M\ corresponding to the CP odd 
scalar. The second matrix gives the masses of the light and heavy Higgs bosons 



= -[ (Ml + m|). ± V(M1 + m|)2 - 4M>| cos^ 2/3 ] . 



The mixing angle a that diagonalizes the matrix (|5.4b| ) can be expressed 

(Ml + m|) 



tan 2a 



rrir, 



tan 2/3 . 



(5.5) 



(5.6) 



If M\ ^ which is the limit where the heavy Higgs is very heavy, this angle coincides 
with /3. This tree level result predicts < Mz- One loop calculations show that this need 
not be the case |2^-^. The third matrix has eigenvalues and Mfj^ corresponding to a 
massless, charged Goldstone boson and a charged scalar. 

Including Eq. ([4.6| ) in the calculations leads to corresponding one loop versions of these 



masses 32,R3l. 



C. Chargino masses 
The following four terms contribute to the chargino masses 



The first two terms are the supersymmetric Yukawa-gauge terms. Letting 
iWi)/y/2, the mass matrix follows 

/ M2 -g2Vd \ / A+ \ 



M2 g2Vu 

V -g2Vd -/i y 



Diagonalization yields two charged Dirac fermions, Ci, C2, with masses 
1 



(5.7) 

iW2 ± 

(5.8) 



Mc^ ,^ = -[(M^ + /i^ + 2m^) T V(^l + + 2m^)' - 4(^2/1 - ml^ sin 2/?)^] . (5.9) 



D. Neutralino masses 



Contributing to the neutralino masses are the terms in Eq. (|5.7| ) and 

- tg'V2^ii+^)KB - tg'V2^li-^)^dB + ^M^BB . 



(5.10) 



The neutralino mass matrix follows 
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-Ml g'vd/V2 -g'v^/^\ ( iB\ 

-M2 -g2Vd/^ g2Vu/^ 

g'vd/V2 -g2Vd/V2 0° 

-g'vjV2 g2Vu/V2 fi ) \ 4>l ) 



(5.11) 



Its eigenvalues are the masses of the four neutrahnos, iVi, iV2, N^, N^. Typically not 
all the eigenvalues of this mass matrix have the same sign. It was noted in that over 
most of the range of the parameters, one of the eigenvalues is of the opposite sign compared 
to the others. This sign is important in distinguishing this neutralino (the "fiippino") and 



especially in deriving mass sum rules El] 



VI. BOUNDARY CONDITIONS AT Mx 



In this paper, as in the previous one [|13[, we work in the modified minimal subtraction 
scheme (MS) of renormalization. The parameters of the Lagrangean are not in general equal 
to any corresponding physical constant. For example, in the case of masses, except for those 
of the bottom and top quark (see [|l3l), all other physical masses, M, will be determined 
from their corresponding running masses by the relation 

M = m(A)|A=M. (6.1) 

This equation is easily solved in the course of an integration of the RGEs for the different 
masses by noting the scale at which it is valid. We have collected the renormalization group 
/3-functions of the MSSM for the gauge and Yukawa couplings to two loops without making 
any approximations in the Yukawa sector in Appendix A. Also included are the one loop 
/3-functions for the soft breaking terms. 

Because SUGRA models make simplifying predictions about the soft parameters at some 
large scale, we initiate the evolution of the renormalization group equations at this scale. 
For simplicity, we take this scale to be the gauge unification scale where we expect the 
gaugino masses to be equal to some common value. It has been demonstrated ||l|] that the 
introduction of supersymmetry leads to gauge coupling unification at approximately ~ 10^^ 
GeV. Therefore we take Mx = 10^^ GeV, and evolve down to 1 GeV, the conventional scale 



at which the running quark masses are given |3^. Furthermore, as already discussed, the 
Higgs potential must be analyzed at some low energy scale that we choose to be Mz- 
At the unification scale, Mx, all the scalars will have a common mass, mo. 



^Q^Mx) = m^^{Mx) = m^,{Mx) = m^^{Mx) 
= m^^{Mx) = m^^Mx) = m^^{Mx) 

= mo , (6.2a) 
the gauginos will also have a common mass, mi/2, 

M,{Mx) = M2iMx) = MsiMx) = m,/^ . (6.3) 
The prefactors of the trilinear soft scalar terms are equal to Aq at Mx 
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TABLE II. Models 

A) mo mi/2 ^0 

General any any any any 

Strict No-scale any 

No-scale any any 

String Inspired any any 

Minimal SUGRA any any any Aq — niQ 



A'^iMx) = A'RMx) = A^i{Mx) ^ . (6.4) 

Also we define the bilinear soft scalar coupling and the mixing mass at Mx by B{Mx) = Bq, 
and iJ,{Mx) = /io- 

Furthermore, to constrain the parameter space, we will take the bottom and tau masses 
equal at Mx 

m,{Mx) = m^iMx) ■ (6.5) 



This being the best motivated mass relation in supersymmetric grand unified theories |^ . 



The purpose of our present analasys is to determine the allowed range of the input 
parameters, Aq, mo, mi/2, and tan (3, which reproduce the known low energy physics. Within 
this general set of solutions, there are four subclasses of soft supersymmetry breaking which 
are of particular interest. Three of these have various soft parameters equal to zero at Mx- 
The fourth predicts a definite relationship among three of these parameters. The first class 
of models follows from the no-scale model and has Aq = mo = Bq = (strict no-scale 



model). In these models, only gaugino masses provide global supersymmetry breaking. The 
second class is the less constraining no-scale case that has only Aq = tuq = 0. A third class 
we consider with Aq = Bq = Q comes from some string derived models. Minimal SUGRA 
models with canonical kinetic terms for the chiral superfields form another class and have 
Bq = Aq — rriQ. Table II lists these four possibilities. In both strict no-scale and string 
inspired cases, we must have -Bo = 0. However, since for us -Bo is an output rather than an 
input variable, -Bq = results must be inferred from its behavior upon varying other input 
parameters. 



VII. NUMERICAL PROCEDURE 

Our numerical procedure relies on the Runge-Kutta technique of numerical integration. 
To make the process of initialization, or the determination of all parameters at one scale, 
more transparent, it is useful to think of the Runge-Kutta integration of the renormalization 
group equations as a vector valued function of a vector variable. The input vector is the 
N-plet of values for all the (N) parameters at Mx- The Runge-Kutta function then returns 
an N-plet result representing the running values of all the parameters at final scales such as 
Mz or 1 GeV. 

In order to make a run (using Runge-Kutta) starting from Mx, we must have the values 
of all the parameters at this scale, but this is difficult to achieve. The problem is that the 
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values of many parameters are known experimentally at low scales. Also, relations ([4.1 7aj) 



and ( |4.17b ) coming from the one loop effective potential hold at such low energies. However, 



the values of other parameters, such as soft breaking terms, are most easily understood at 
higher energies where theoretical simplification {e.g., universality) may be invoked. We 
therefore have no scale at which there is both theoretical simplicity and experimental data. 

We can phrase the problem in another way. Choosing Mx (where we have theoretical 
simplicity) as our starting point, can we find the Mx values of all parameters such that we 
recover the expected low energy values after renormalization group evolution to experimental 
scales, and all constraint relations among the parameters are satisfied at the appropriate 
scales? We distinguish between constrained and free parameters. The former are constained 
by experiment {e.g., quark masses, gauge couplings, etc.) or relations among themselves 
{e.g., bottom r Yukawa equality at Mx, minimization conditions at Mz, etc.). Given the 
present experimental data, the latter cannot be constrained by these two criteria and must 
be viewed as input parameters. After we have run, we find some ranges of their values to be 
inconsistent with phenomenological considerations such as electroweak symmetry breaking 
and the resulting sparticle spectrum. The exact values of the constrained parameters at Mx 
are affected by the choice of values for the free parameters at Mx since the evolution of all 
parameters are coupled. Therefore, given a choice for the free parameters at Mx, we must 
find the Mx values of the constrained parameters consistent with all the constraints. 

The functional nature of the (Runge-Kutta) integration allows us to define a set of n 
(where n is the total number of constrained parameters) equations 

Gfc(xo) = , (7.1) 

where k = l,...,n, and xq represents the values of all constrained parameters at Mx- 
The solution to the initialization problem, therefore is reduced to solving a system of n 
simultaneous non-linear equations in n unknowns. In this work as in Ref. |T^, we will use 
routines, based on the "shooting" method, to solve systems of nonlinear equations. The 
method involves making a guess for the solution, then assessing its merits based on how 
well the equations are satisfied, given some tolerance. The next guess (or shot) is adjusted 
according to how accurate the previous one was. The process is optimized and iterated until 
the routine converges on a solution. In this way, the values of all parameters are ascertained 
at one common scale which we take to be Mx- When the parameters are evolved to lower 
scales using these initial values, they attain the experimentally known values and they satisfy 
any relations amongst themselves that were used as constraints in the shooting procedure. 

As discussed previously, we start our runs at Mx at which scale we can make simplifying 
assumptions about the soft breaking terms based on various SUGRA models. This requires 
that we use the solution routines to consistently find the Mx values of all known low energy 
parameters such as lepton and quark masses and mixing angles and gauge couplings. This 
amounts to solving for sixteen unknowns (nine masses, three angles and a phase, and three 
gauge couphngs). Alternatively, we could start our runs at Mz or 1 GeV; however, this 
now requires solving for sixty-three unknowns (the values of the soft breaking terms at low 
energy) that must evolve to just four different values at Mx. The efficiency of the former 
method is obvious. 

There are ostensibly seven parameters in the models we consider. These are Aq, Bo, mo, 
^^1/2, f^o, tan/3, and rrit. The two minimization constraints ( |4.13| ) and ( [4. 14] ) reduce this 
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set to five, which are taken to be Aq, mo, mi/2, tan/3, and rrit. In the present framework, 
Bq and /io will be determined using the numerical solutions routines in conjunction with 
the minimization of the one loop effective potential at Mz in the process of evolving from 
Mx to 1 GeV. Minimization at Mz will give B{Mz) and ii{Mz)- To arrive at Bq and Hq 
(their corresponding values at Mx), we employ the solution routine as follows. A guess 
for Bq and fiQ is made at Mx and then the parameters of the model are run to Mz at 
which scale the evolved value of B is compared to the minimization output value for B at 
Mz- The same is done for /x. If the compared values agree to some set accuracy, then Bq 
and fiQ are the required values. Other analyses that also extract B{Mz) and fJ,{Mz) simply 
evolve these two parameters via their renormalization group equations back to Mx to find 
Bq and /io relying on their near decoupling from the full set of RGEs. We note that the 
sign of /i is not determined from the minimization procedure, thus we must make a choice 
for it. To constrain the parameter space further, the bottom quark and tau lepton masses 
will be taken equal at Mx- This equality is a characteristic of many SUSY-GUTs. This 
constrains the model to four free parameters, Aq, mo, mi/2, and tan/?. Demanding that 
mb{Mx) = mr{Mx) and achieving the correct physical masses for the bottom quark and 
tau lepton fixes the mass of the top quark which affects the evolution of the bottom Yukawa 
significantly. We shall assume gauge coupling unification, an assumption which appears 
reasonable when one considers SUSY models with SUSY breaking scales ^ 10 TeV. 

In a complete treatment, the solution routines would be used to find the precise (similar) 
values of ai, 02, and 03 at Mx that will evolve to the experimentally known values at Mz, 
however this increases the CPU time considerably. We shall therefore sacrifice some precision 
in their Mz values by taking them exactly equal at Mx- This is already a theoretical 
oversimplification since one does not expect the gauge couphngs to be exactly equal due 



to threshold effects at the GUT scale We find that for all cases we have studied, the 
common value a{^{Mx) = a2^{Mx) = a^^{Mx) = 25.31 leads to errors no bigger than 1%, 
5%, and 10% in ai{Mz), a2{Mz), and a^^Mz), respectively. This is not so bad considering 
that the (combined experimental and theoretical) errors on a^lMz) from some processes 
can be as large as 10% 

It is well known that there is a fine tuning problem inherent in the radiatively induced 
electroweak models. For certain values of the parameters, the top quark mass must be 
tuned to an "unnaturally" high degree of accuracy to achieve the correct value of Mz- This 
problem is generally handled by rejecting models that require "too much" tuning. The 
amount of tuning is usually defined quite arbitrarily. The usual procedure is to define fine 
tuning parameters 



x2 9M| 



(7.2) 



where Xi are parameters of the theory such as mo, m.1/2, Z^, or rrit- One then demands that 
the Ci be less than some chosen value that is typically taken to be 10. 

We have analyzed to some extent the differences in using the tree level vs- one loop 
effective potential. The basis for the "theoretical" fine tuning problem can be seen (for 
example, assuming the renormalization group equation for the up sector Higgs mass is 
dominated by the top quark Yukawa coupling) in the dependence of M^y on the top quark 
Yukawa coupling yt |^ 
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(7.3) 



We remark that this fine tuning problem is exacerbated if one uses only the tree level analysis 
of the potential. The vacuum expectation value coming from the minimization conditions 
of the tree level potential changes rapidly from to infinity over the interval (A^, A^,). Using 
the prescription of Ref. ||28| for the scale A at which to adequately minimize the tree level 
potential, to extract f (A), and thereby to arrive at a value for M^, one finds that although 
a small variation in yt{Mx) may lead to a small variation in A;,, the steepness in the tree 
level VEV can lead to a large variation in the value of v{A) and therefore in Mz- Hence, in 
the tree level analysis, solutions which may be within the bounds of the "theoretical" fine 
tuning may nevertheless display a fine tuning aspect because of this "tree level" fine tuning 
of A. 

However, our use of the one loop effective potential (i.e., including AVi) stabilizes the 
VEV around the Mz scale and this particular fine tuning goes away. The true VEVs depend 
on scale through wave function renormalization effects which are never large as can be seen 
from the form of the renormalization group equations for the VEVs in the Appendix A. 

In this analysis, we shall also reject solutions based on fine tuning considerations; how- 
ever, our method differs somewhat from the usual one in that it is incorporated in the 
solution routine described above. The routine is an iterative one which determines the con- 
vergence properties of the solution. Very slow convergence reflects an inherent flne tuning. 
Therefore, if the convergence is too slow, we will reject the solution. Effectively we are 
rejecting any solution which the computer cannot pinpoint within an allotted number of 
iterations. 

Given values for Aq, rriQ, mi/2, tan/3, and sign(yu), the solution routines search for the 
values of v{Mx), mf,^r{Mx), rrit^Mx), Bq, and /xq- The process by which Bq and fiQ are 
found was described above. The remaining three parameters are determined similarly. The 
routine makes a guess for v{Mx), mb^r{Mx), and mt{Mx), then the full renormalization 
group equations are evolved to 1 GeV calculating superparticle threshold masses in the 
process and minimizing the one loop effective potential at Mz- The merits of the guess 
for v{Mx), TTib^riMx), and mt{Mx) are assessed by comparing the resulting values of M^, 
rriril GeV), and mi,{l GeV) with the expected ones. The process is iterated until the correct 
values are achieved to within a tolerance of 1%. 

By a "solution," we will mean that a choice for the four free parameters (plus a choice for 
the sign of /x) at Mx results in a complete set of all parameters at Mx (using the shooting 
routines) consistent with electroweak breaking {Mz = 91.17 GeV), with equality of the 
bottom and r masses at Mx, and with low energy experiment. 

Such a solution then yields a precise spectrum of sparticle masses. This spectrum can 
be further used to limit the allowable free parameter space by subjecting the spectrum to 
experimental restrictions. 

VIII. THRESHOLDS 

In the minimal low energy supergravity model being considered, the super particle spec- 
trum is no longer degenerate as in the simple global supersymmetry model in which all the 
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super particles are given a common mass, Msusy- In this simple case, one makes one course 
correction in the renormalization group evolution at Msusy- In the model with soft symme- 
try breaking, however, the nondegenerate spectrum should lead to various course corrections 
each occuring at the super particle mass thresholds. To this end, the renormalization group 
/^-functions must be cast in a new form which makes the implementation of the thresh- 
olds effects (albeit naive) evident (see Appendix B). Since the MS renormalization group 
equations are mass independent, particle thresholds must be handled using the decoupling 



theorem |^T|, and each super particle mass has associated with it a boundary between two 
effective theories. Above a particular mass threshold the associated particle is present in 
the effective theory, below the threshold the particle is absent. 

The simplest way to incorporate this is to (naively) treat the thresholds as steps in the 
particle content of the renormalization group /^-functions This method is not always 

entirely adequate. For example, in the case of the SU (2) gauge coupling there will be scales 
in the integration process at which there are effectively a half integer number of doublets 
using this method. For example, in the region Mb < A < Mt between the mass of the 
bottom and top quarks, the number of quark doublets is taken to be two and a half using 
this method. A similar situation occurs with sparticle doublets. We believe, nevertheless, 
that this method does yield the correct, general behavior of the evolution. It is a simple 
means of implementing the smearing effects of the non-degenerate super particle spectrum. 
The determination of the spectrum of masses is done without iteration as is common in 
other analyses. Our method deduces the physical masses by solving the equation m(A) = A 
for each superparticle in the process of evolving from Mx to 1 GeV. The usual iterative 
method requires several runs to find a consistent solution. 



IX. ANALYSIS 

The tremendous computing task involved in analyzing the full parameter space of the 
soft symmetry breaking models, using the methods described as designed, would be far too 
time consuming given the computing facilities available to us. Therefore, in the following 
analysis, some simplifications will be made in the procedural method. First, only the heaviest 
family of quarks and leptons will have non-zero mass. This will cut down on the CPU time 
required for the solution routines to consistently find the running mass values at Mx since 
there are three instead of nine masses to determine after this simplification. Second, as 
stated previously, the value of the strong coupling at Mz will be allowed to vary from its 
central value of .113 by at most 10%. This translates into a similar error in the bottom quark 
mass. Third, the allotted number of Runge-Kutta steps, involved in numerically integrating 
the renormalization group equations, will be cut down to ~ 100. 

Our method involves four input parameters Aq, mo, mi/2, tan/? (and the sign(/i)). The 
output is Bq, /io, Mt, M/j, and all the masses of the extra particles associated with the 
MSSM. Efficient use of CPU time required that we proceed as follows. For a given model, 
our initial exploration of the parameter space was performed in a coarse grained fashion, 
tan/? was most commonly coarse grained as 2, 5, and 10, with only some runs involving 
higher values (e.g., 15, 20). The other three input parameters were varied in steps of 50 and 
100 GeV. Values larger than ~ 500 GeV were rarely ever used. We subsequently narrowed 
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down on the allowed hyperenvelope by fine graining around the edges of the expected region 
(based on the coarse graining results). 

Our raw data consists of those runs which satisfy the following two criteria. The first 
is consistency with electroweak breaking; that is, the correct value of Mz is achieved from 
the minimization of the one loop effective potential with fi'^{Mz) > 0. The second criterion 
is no fine tuning, as implemented in our method (see Section VII). The solution routines 
employed return a numbered code representing the convergence properties of the solution 
which we use to screen the runs. 

The raw data can then be progressively filtered based on at least three physical con- 
straints. The first is cosmological. Conserved R-parity (+1 for particles and —1 for super- 
particles) requires the existence of a stable lightest supersymmetric particle (LSP). Astro- 
physical considerations indicate that the LSP must be neutral and colorless. Cosmological 
considerations based on the LSPs contribution to the density of the universe indicate that it 
must have a mass less than ~ 200 GeV [^3]-^. Therefore, points in parameter space that 
lead to LSPs other than neutralinos with masses less than ~ 200 GeV are cut. Second, flavor 
changing neutral current bounds can be used to reject runs in which interfamily splitting 
of squark and slepton masses is too large, although this is seldom the case in the type of 
low energy supergravity models with soft breaking universality that we consider. Third, 
experimental limits on the masses of the superparticles can also be used as another criterion 
to reject some runs. 

X. RESULTS 

In the following, we discuss the results of our analysis of the four classes of models listed 
in Table II and of a general, non-constrained class. In the general case, we discuss results 
based on all our runs with no constraints on the soft parameters. Because we have performed 
a coarse grained study, using fine graining only for selected regions, our results should only 
be considered quahtatively valid. Based on these, we hope to be able to ascertain the general 
trends in the data, and make some general predictions about the feasibility of the models 
considered. 

Because the GUT inspired constraint Eq. (|6.5| ) is enforced in this analysis, the results 
will depend on the mass of the bottom quark. Namely, lower bottom quark masses require 
larger values of the top quark mass to satisfy this relation, and higher bottom quark masses 
require smaller values of the top quark mass. For simplicity, most results will be reported 
for the case mfe(l GeV) = 6.00 GeV, but lower mass (5.70 GeV) and higher mass (6.33 GeV) 
cases were also studied. These lower and higher values are within experimental uncertainty. 
The running value of 6.00 GeV for mb(l GeV) corresponds to a physical bottom mass 
Mb = 4.85 ± .15 GeV, with the uncertainty coming from the error in the strong coupling, as 
discussed above. 

In Table III, we present the allowable ranges of the soft breaking parameters for the 
general unrestricted case and for each of the four subcases considered. Since Bq and fio 
are not input parameters in our method, we must infer the results for the strict no-scale, 
string inspired, and minimal SUGRA cases. We do this by setting a small tolerance on |i?o| 
and \x\ of 25 GeV. In compiling this table, we have discarded all runs with Mf < 120 GeV 
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TABLE III. Allowed Soft Parameter Ranges (in units of GeV). 



Case Aq mo mx/2 -Bp 

General (0,500) (0,800) (190,400) (33,817) (-56,1012) 

No-Scale (190,262) (40,230) (-60,130) 

Strict No-Scale (210,240) (159,188) l^ol < 25 

String Inspired (0,200) (200,250) (88,188) i^oi < 25 

Minimal SUGRA (0,500) (0,300) (200,300) (159,467) (-12,407) 



based on what we feel is a reasonable experimental lower bound for the top quark mass. We 
note from the table that /io is less than ~ 200 GeV in the first three cases, and only in the 
minimal SUGRA case does it reach larger values. The largest value of /iq was 817 GeV and 
occured for a run that does not correspond to any of the above cases. The values of the 
other soft parameters in this run and the resulting spectrum are displayed in Table IV. 



A. General Case 



Based on all runs and imposing no particular constraints on any soft parameters, we find 
that 86 < < 185 GeV for the top quark and 35 < Mh < 141 GeV for Higgs boson. Table 
IV contains three scenarios with input parameters that do not conform to any of the four 
aforementioned cases. Scenario (a) has a very large value of rriQ. The sparticle masses are the 
largest we have encounterd so far in our discussion. One interesting feature of this scenario 
is that the lighter stop is lighter than any of the sleptons. Scenario (b) has a gluino of mass 
908 GeV and a top quark of mass 184 GeV. The Higgs boson mass is one of the heaviest of 
all the scenarios discussed at 132 GeV. Both (a) and (b) have very large neutralino masses 
~ 170 GeV. The spectrum in (c) represents another scenario in the general case, but with a 
larger value of tan /5. 



B. No-scale Case 



In all cases considered, the mass of the LSP, when it is a neutralino, is observed to be 
correlated with the value of mi/2- Therefore, we find that mi/2 cannot be taken too large 
(^ 400 GeV). In the no-scale case, we present plots of Mt vs. mi/2 for three values of 
mb{l GeV) and containing all points satisfying the various criteria outlined in Section IX. 
Figure 1 is such a plot for mb{l GeV) = 6.00 GeV. The "right edge" of the envelope is 
defined by points whose neutralino LSPs are just slightly heavier than the lightest charged 
superparticle (usually a for us). The "left edge" defines the threshold of consistent 
electroweak breaking. The top and bottom edges are set by the requirement that Eq. ( |6.5| ) 
hold. The range of tan/3 considered leads to definite lower and upper bounds on Mf p|-[l2|. 
In Figs. 2 and 3, we display similar plots with lower (mb(l GeV) = 5.70 GeV) and higher 
(mfe(l GeV) = 6.33 GeV) bottom quark masses. The previously noted dependence of the 
data on Mf, is evident from these figures. 
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TABLE IV. Sample Particle Spectra for General Models (in units of GeV) . 

(a) (b) (c) 

^0 200 

mo 800 500 100 



400 400 300 



tan/3 2 5 15 

sign(/Lt) + + + 

Bo 638 326 44 

/xo 817 687 301 

Mt 175 184 155 

g 913 913 692 

ui,U2 1112,1137 926,956 604,629 

tiM 696,993 626,869 476,636 

Ji.da 1109,1139 922,959 602,634 

61,62 966,1109 831,923 573,611 

61,62 815,847 525,573 159,236 

fi,f2 815,847 521,573 143,241 

844 568 223 

Ur 844 567 221 

Ci 324 323 222 

C2 724 604 358 

TVi 171 171 125 

N2 324 324 223 

iVg -712 -590 -326 

Ni 726 604 358 

h 118 132 112 

H 1149 755 178 

A 1147 755 177 

1150 759 194 
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From these, we note that the available range of mi/2 decreases with increasing Mi,. Fig. 1 
indicates that 190 GeV ^ "021/2 ~ 265 GeV. We can draw no conclusions, however, about 
the value of tan (3. Although, we found that values of tan (3 larger than 18 never led to 
solutions in the no-scale case. We can conclude from the mb{l GeV) = 6.00 GeV case that 
85 ^ Mt < 132 GeV. Prom similar plots involving Mh, we conclude that 35 ^ Mh ^110 
GeV. In Fig. 4, we display the familiar dependence of Mf on tan /3 for a particular value of 
mi/2 in the allowed range. 

The official experimental lower bound on the top quark mass is 108 GeV. The figures 
indicate that the top quark cannot have a mass greater than ~ 132 GeV in this model, if 
mb{l GeV) = 6.00 GeV. This upper bound is raised to ~ 160 GeV, if mb{l GeV) = 5.70 
GeV, and the model is ruled out, if mb{l GeV) = 6.33 GeV. 

Figure 5 is similar to Fig. 1, but we have chosen the sign of fi negative. The allowed 
region is displaced down with respect to the positive 11 case with the upper bound on 
now 113 GeV and very close to the experimental limit. 

Table V displays various spectra of superparticle masses in the no-scale case. In spectrum 
(a), we present a scenario in which mi/2 = 254 GeV and tan/3 = 5. In this particular 
scenario, the top quark mass is 131 GeV and slightly below is the mass of the light Higgs 
boson at 98 GeV. A qualitative feature of most spectra is that the sleptons are lighter than 
the squarks. Another feature is that one of the stops is lighter than all other squarks. In 
this table, as in similar ones to follow, we do not include the second family sfermion masses 
since they are generally degenerate with those of the first family. Also, we have singled out 
the flippino by associating with it a negative mass. The input parameters of spectra (b) 
and (c) difi'cr only in the respective values of tan f3, 3 and 10. The spectra are identical in 
almost every respect except for the splitting in the stop masses. The stop mass splitting in 
spectrum (b) is 154 GeV and 140 GeV in spectrum (c) as expected since the product //cot P 
is bigger in case (b). Spectrum (a) represents a shifted version of (b) or (c). 

C. Strict No-scale Case 

The results of the strict no-scale case must be interpolated from the no-scale results, 
because Bq is not an input parameter in our procedural method. Therefore, we plot in Fig. 6 
Mt vs. Bo and deduce the Mt bounds from slicing the data along Bq — 0. Admittedly, the 
relatively small number of points makes the perimeter of the region unclear in some areas. 
We get approximately 90 < Mt < 127 GeV. For M,„ we find 57 < M;, < 110 GeV with an 
uncertainty in the lower bound of ~ 10 GeV due to lack of definition in the lower end of 
the envelope. Inspection of the data indicates that 3.5 ^ tan/3 ^ 9. Thus, it appears that 
tan P cannot be too small or too large to accommodate the strict no-scale case. This case is 
a special case of both the no-scale and string inspired cases. 

Finally, in Table VI, we display spectra of superparticle masses for three representative 
strict no-scale scenarios with various tan/5 and mi/2- Although, as evident from the table, 
none of these scenarios have Bq strictly equal to zero, they do satisfy |i?o| < 25 GeV, which 
is the tolerance we have settled upon. The LSPs are neutralinos in these three scenarios 
and have masses that increase with increasing mi/2. The top quark masses are comparable 
at ~ 126 GeV and just above the experimental hmit, and the Higgs bosons masses are also 
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TABLE V. Sample Particle Spectra for No-Scale Models (in units of GeV). 

(a) (b) (c) 

^0 

mo 

mi/2 254 240 240 

tan/3 5 3 10 

sign(//) + + + 

Bo 29 66 -12 

Ho 203 198 166 

Mt 131 126 124 

g 589 558 558 

ui,U2 508,530 482,502 482,501 

ti,t2 408,559 381,535 392,532 

di,d2 507,535 481,507 481,507 

61,62 507,508 476,486 472,491 

61,62 107,183 101,173 102,174 

fi,f2 105,184 100,173 97,176 

i>e 166 157 155 

i>r 166 157 155 

Ci 164 150 145 

C2 296 293 267 

iVi 100 90 92 

N2 170 157 152 

N3 -246 -240 -208 

N4 298 296 268 

h 98 88 99 

H 278 290 203 

A 276 285 202 

287 296 217 
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approximately equal ~ 99 GeV. As in the no-scalc case, we observe that increasing the value 
of mi/2 has the effect of shifting the squark and slepton spectra up. 

D. String Inspired Case 

As in the strict no-scale, our results for the string inspired case necessitates interpolating 
^0 = data to Bq — 0. Once again we find the perimeter of the allowed region is not well 
defined everywhere. Hence, our results are only qualitative. Figure 7 is similar to Fig. 6, but 
in this case we only fix Aq = 0. Since the strict no-scale case is a special case of the string 
inspired one, the 127 GeV top quark mass upper bound is not expected to decrease but 
rather to increase in this case. Shcing along Bq = yields 85 < < 140 GeV. Similarly, 
for the hght Higgs we get 57 ^ Mh < 113 GeV. The data indicates in this case, as in the 
strict no-scale case, that there is a lower bound on tan /3 of ~ 3. 

In Table VII, we collect some spectra in the string inspired case. As in the strict no-scale 
case, we have settled on a tolerance of |-Bo| < 25 GeV for the scenarios considered. Spectrum 
(a) with input parameters Aq = 0, mo = 100 GeV, mi/2 = 200 GeV, and tan/? = 10 has 
Bq — —4 GeV. The strict Bq — run has a value of tan/3 between 5 and 10, but the 
qualitativcc features should be very similar to the scenario we present with tan/5 = 10. 
Scenario (b) is similar to scenario (a) except it has a larger value of mi/2 = 250 GeV. As 
noted before, this leads to an upward shift in the squark and slepton spectrum of (b). The 
upward shift is 23% in the squark masses and 15% in the slepton masses. Scenario (c) is 
similar to scenario (a) except it has a larger value of mo = 200 GeV. This shifts up the squark 
and slepton masses of (c) with respect to (a). In this case, the squarks shift up by only 7% 
and the sleptons shift up by 49%. Therefore, spectra (b) and (c) have the curious property 
that the squarks of (b) are heavier than those of (c), but the sleptons of (b) arc lighter than 
those of (c). Lastly, we observe that the lighter selectron and stau of all spectra discussed so 
far display a comparatively large mass gap from the rest of their respective slepton partners, 
however in scenario (c) this mass gap is comparatively small; and all sleptons are almost 
degenerate. 

E. Minimal SUGRA Case 

We now consider models for which the relation Bq = Aq — mo among the soft couplings 
holds. We find it convenient to define x = ~ Bq — mo, which should equal zero in this 
case. In Fig. 8, we plot Mt vs. x- Taking a slice along x = of the region depicted in this 
plot, we find that 82 ^ Mt ^185 GcV. A similar plot was used to arrive at bounds for the 
light Higgs boson 57 ^ Mh ^ 139 GeV. Our analysis indicates that there is no preferred 
values of tan /? in this case. 

Table VIII includes the spectra for some particular scenarios in the minimal SUGRA 
case. The three spectra appearing in this table have the interesting feature that the slepton 
masses are larger than they have been in all cases discussed so far. In fact the slepton masses 
have values comparable to the lighter stop in (a), (b), and (c). These three scenarios each 
have relatively large, non-zero values of Aq. As expected, scenario (b) has the larger value of 
II cot P and therefore has the larger splitting in stop masses when compared to (a). Indeed, 
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TABLE VI. Sample Particle Spectra for Strict No-Scale Models (in units of GeV). 



(a) 



(b) 



(c) 



^0 

mo 

■mi/2 

tan/3 

sign(/i) 

Bo 

Mo 






240 
8.3 

+ 
-4 
171 





230 
10 

+ 
-11 
159 






210 

5 

+ 

18 
159 



Mt 
9 



126 
558 



124 
536 



128 
491 



Ui,U2 
tl,t2 

di,d2 
bi,b2 

61,62 
fl,f2 



482, 501 
390, 533 
481, 507 
473, 490 
102, 174 
98, 175 

156 

155 



462, 481 
373, 514 
462,487 
453, 472 
99, 167 
93, 169 

148 

148 



424, 441 
331, 480 
423, 447 
418,430 
92, 154 
90, 154 

133 

133 



C2 



147 
270 



137 
259 



122 
253 



iVi 

iV4 



92 
153 
-213 
271 



87 

144 
-200 
260 



76 

131 
-197 
255 



h 
H 

A 



99 
222 
221 
235 



99 
194 
193 
209 



96 
223 
221 
234 
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TABLE VII. Sample Particle Spectra for String Inspired Models (in units of GeV). 

(a) (b) (c) 

^0 

mo 100 100 200 
■mi/2 



200 250 200 



tan/3 10 10 10 

sign(/Lt) + + + 

Bo -4 -6 20 

lio 132 168 88 

Mt 124 125 123 

g 469 580 469 

ui,U2 416,431 510,530 449,464 

ii,i2 330,467 417,556 361,489 

di,d2 416,438 509,536 449,470 

61,62 407,424 509,509 440,453 

61,62 134,178 145,206 219,248 

fi,f2 129,179 141,208 216,248 

i>e 160 191 235 

i>r 159 190 234 

Ci 109 150 77 

C2 231 272 214 

Ni 71 96 53 

^^^2 119 157 105 

N3 -170 -209 -120 

N4 232 273 214 

h 98 100 98 

H 189 226 230 

A 188 225 229 

H"^ 204 238 242 
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TABLE VIII. Sample Particle Spectra for Minimal SUGRA Models (in units of GeV). 

(a) (b) (c) 

^0 400 500 400 

mo 200 200 300 

■mi/2 



300 300 200 



tan/3 10 5 10 

sign(/Lt) + + + 

Bo 212 288 115 

lio 404 436 158 

Mt 164 163 124 

g 692 692 469 

ui,U2 627,652 627,652 500,513 

ii,i2 415,643 389,644 330,529 

di,d2 625,656 625,656 500,519 

61,62 573,630 574,627 463,502 

61,62 235,293 235,292 313,334 

fi,f2 221,295 230,294 306,333 

i>e 282 282 324 

i>r 279 281 322 

Ci 236 236 123 

C2 438 471 247 

iVi 127 126 76 

N2 236 236 129 

N3 -421 -454 -200 

N4 438 472 248 

h 119 118 101 

H 371 491 325 

A 371 490 325 

H"^ 379 497 334 
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this is the only significant difference in these two scenarios. Scenario (c) in the minimal 
SUGRA case has a fairly degenerate slepton sector like that of (c) in Table VII. Once again 
we observe the curious feature when comparing (a) and (c) that the squarks of (a) are 
heavier than those of (c), but the sleptons of (a) are lighter than those of (c). It appears 
that the run with the larger value of mo has the heavier sleptons and that with the larger 
value of mi/2 has the heavier squarks. Scenarios (a) and (b) are interesting because they 
have large values of /xtan/? which leads to large splittings in the bottom squark and stau 
masses. This represents an example where neglecting the bottom or r Yukawa can radically 
alter the spectrum. The value of Mf in (a) and (b) is ~ 163 GeV, however in (c) it is 124 
GeV. We find based on studying other similar runs that the input parameter accounting 
mostly for this large difference is the value of mi/2. Both (a) and (b) have mi/2 = 300 GeV 
whereas (c) has mi/2 = 200 GeV. We observe also that (a) and (b) have larger (comparable) 
values of /iq and Bq whereas (c) has a ~ 50% smaller values of these input parameters. 



XI. CONCLUSIONS 

Minimal low energy supergravity models were considered. They have the appealing 
feature that the electro-weak symmetry is radiatively broken for certain ranges of the soft 
breaking parameters and of the top quark mass. The study of specific models, with some soft 
parameter fixed or related, resulted in upper bounds for the top quark mass. No-scale models 
in which only gaugino masses provide global supersymmetry breaking yield top quarks with 
masses less than ~ 127 GeV. The results are sensitive to the value of the bottom quark 
mass. Lower bottom quark masses, within the experimental uncertainty, lead to higher top 
quark upper bounds. In these models, the ratio of vacuum expectation values of the two 
Higgs fields is expected to be larger than ~ 70°. 

Although the perimeter of the allowed regions were often fuzzy, we could, nevertheless, 
draw some general conclusions from our results. For all our runs, with no restrictions on the 
soft terms, we find for the top quark Mt ^ 185 GeV and for the light Higgs boson ^ 141 
GeV . 
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APPENDIX A: THE SUSY /3-FUNCTIONS 



Using some of the notation of Falck the superpotential and soft symmetry breaking 
potential are as follows: 
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W = tY„$„g + dY^^^Q + #Y,|./ + + h.c. , (Al) 



3 



V^ausmo ^ ^ E ^AAi + ^-C. . (A3) 
^ Z=l 

Various (J2's have been omitted and a sum over the number of famihes is imphed in the 
squark and slepton mass terms. Also, hats imply superfields and tildes the superpartners of 
the given fields. 

We start with the gauge couplings 

II ^ "16^^'^' + (1^ % ~ Tr{Q„Yt Y„ + QdYlY, + QeYtY j] (A4) 

where i = InA and I — 1, 2, 3, corresponding to gauge group SU(3)c x SU(2)l x U(1)y of 
the Standard Model. 

In the Yukawa sector the /3-functions are 

The evolution of the vacuum expectation values of the Higgs's is given by 



1. The one loop SUSY /3-functions 

The various one loop coefficients for the gauge couplings are defined to be 

' bi = -|-2ng , 

< 62 = 5-2ng , (A7) 
^ 63 = 9 - 2n, , 

with Ug = \nfi. 

In the following, we list the one loop contributions for the parameters of the superpo- 
tential. 

^ = 16^t Tr{3YtY„ + 3YtY, + Ytyj - Z{\gl + gl) ] . (A8) 
The one-loop contributions for the Yukwas are given by 
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f^i'^ = 3YtY„ + YtY, + 3Tr{Yt YJ - {^^gf + 3gl + fgl) , (A9) 

/?(^) = 3YtY, + Yt Y„ + IV{3YtY, + YIyJ - {^gf + Sg^ + ^gl) , (AlO) 

/?(^) = SYtY, + Tr{3YtY, + Y^YJ - {^gl + Sgj) . (All) 
The one-loop contributions for the running VEVs are 

7^ - ^(^^1 + 92) - 3Tr{Yt YJ , (A12) 

72 = l^lal + gl) - 3Tr{YtY,} - Tr{YtYj . (A13) 
The soft symmetry breaking terms are known to us only to one loop 



fjAij 1 ykj yik Aij 

-^-Y6^2i ^{^e^rA^^^ + 5Af ^(YtYj'^^- - 3^(Y,YtYj^^- 

+ 2{Al"'\Y^"'\^ + 2,A'r\Yt'f) - QilglM, + ^^Ma) ] , (A14) 

5 

yik ykj 

+ (4' - 4')(YlYJ^^-^ + 2(Y,Yt)^'=A^^-^ + 2{A^-\Y^-\' 

14 S2 
+ 3^r|yr n - ^^^i^i - 6^2'M2 - -glMs ] , (A15) 

^ = 4(Y„Yt)^^yl^^-^ + 5yi;^^(YtYj'=^- - 3^(Y„YtYj'^- 

yik ykj 

+ (^^ - ^^^KY^Y,)^^-^ + 2{Y^YlfA'/^ + 6At-|rr r 

26 32 

- -glM, - 6glM, - -glM^ ] , (A16) 

^ = ^[ E 3|i;^f + + m^^. + jA^f ) + ^^^TrlFm^ - ^gM 

ij 

- ^glMl ] , (A17) 
^ = ^[ E(ln^f « + < + + l^if ) + m'Wrnl^ + ml + ml 

ij 

+ l^rn) - ^^iTr{ - ^gfMl - SglM^ ] , (A18) 

^ = ^[ E2|n^^r« + < + < + l^^f) + ^^?Tr{rmn - fgM ] , (A19) 

^ = ^[ E in^T « + < + < + l^if) - ^^i^Tr{Fmn 
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- ^gfM', - ?>glMl ] , (A20) 

- ^/iM! - fglMl ] , (A21) 

- ^^gfM! - fglMl ] , (A22) 



dt Stt , , 



+ m')) + ^gfTriYm^} - ^gfM^ - SglM^ - fglM^ ] , (A23) 

^ = ^[ SA^^irj^f + SA^\Y^^ + A^\Y:^^ - ^gfM^ - Sg^M, ] . (A24) 

where, as in Falck, sums are implied over all indices not appearing on the left hand side and 
where 

ng 

TT{Ym^} = J2i^Q, - 2"^i + ml - m\ + m^J + - . (A25) 

i=l 

Note that in an anomaly free theory, this term is zero if all the masses are equal at some scale. 
That it is zero at such a scale can be seen by using the fact that there is no gravitational 
anomaly. In such a case, Tr{ym^} = m^Tr{y} = 0. To show that it remains true at 
all scales, one also needs the cancellation of the other triangle anomalies. For example, in 
Eqs. ( |A17| - |A23|) , the g\M\ terms come from one loop mass corrections involving two bino- 
particle-sparticle vertices and are therefore proportional to . Thus one needs to have the 
f/(l) anomaly cancellation condition Tr{y^} = in order to show dYx\Yvp?^jdt = 0. The 
gaugino masses evolve as follows 

din Ml 1,9 s 
^ = ob,gf ■ A26 



2. The two loop SUSY /3-functions 



The two loop contributions to the gauge couplings are 



'Ik 





38 


6 


88 \ 




15 


5 


15 






2 

5 


14 


8 




V 


11 

15 


3 


68 
3 / 





_9_ 

25 
3 

5 






(A27) 



and 
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26 14 18 

5 5 5 

(Cif) = I 6 6 2 I , with / = u ,d ,e , (A28) 
4 4 

The two-loop contributions to the Yukawa couphngs are given by 

= -4(YtYj2 - 2{YlY,)' - 2YtY,YtY„ - 9Tr{YtYjYtY„ 

- Tr{3YiY, + YtYjY^Y, - 3Tr{3(Yt YJ^ + Y^Y^Y^YJ 

+ ilgl + QglM^u + ilgD^l^, + i^gl + ^Qgl^n^u} 

,26 403, 4 , 21 4 ,32 304, , 

+ ghl + ^gfgi + ^ghi (A29) 

15 

/3f = -4(YtY,)2 - 2(Yt YJ2 - 2YtY,YtY, - 3Tr{YtYjYt Y„ 

- 3Tr{3YtY, + YtY^Y^Y, - 3Tr{3(YiY,)2 + (YtYj^ + Y^Y^Y^ YJ 

+ ilglM^^ + i^g! + QgD^lr^, + (-^^? + i6^3^)Tr{YtY,} + (^^;?)^{y1yj 

+ ^h""' + 18^^' ^ ^^""^ ~ T^^^' ^ " + + 1^'^' ^ ^^'^^'^ ^^^^^ 

/3f ) = -4(YtYj2 - 3Tr{3YtY, + YlYjYtY, - 3Tr{3(YlY,)2 + (YtYj^ 

+ YtY.Yt YJ + i6gl)YlY, + (^^??)Tr{YtYj + {-^gf + 16^?32)Tr{YiY,} 

+ (f % + ^)gt + (6% - y )^72 + Ighl (A31) 
The two-loop contributions to the anomalous dimension of the scalars are 

7g = ^Tr{3(Yt YJ^ + 3Yt Y^yIy,} - {^gl + \gl + 20^7DTr{Yt YJ 

,279 1803 , 4 ,207 357 , 4 ,27 9 ,3 . ^ 

- W + 1600"^)^^ - + "6?"^)^^ - ^80 + 80"^)^^^^ ' ^^'2^ 

7g = ^Tr{3(YtY,)2 + 3YtY,Yt Y„ + (YIyJ^} - {\gl + ^gl + 20^32)Tr{YtYj 

,9 2 3 r^^t^^ T ,279 1803 ,4 ,207 357 , 4 

- -9i + -99)Tr|YJY 1 - ( \ - ( \ njg^ 

igyi^2^2; \ e eJ ^800 ^ 1600 ^^^^ ^ 32 64 ^'^^ 

27 9 

^80 80 ^'^^^^ ^ ' 



APPENDIX B: THRESHOLDS 



Ref. 1^ gives formulas valid to two loops to compute the /3-functions of gauge, Yukawa, 
and scalar self-quartic couplings in a general gauge theory. These will be useful in obtaining 
the required form of these /3-functions for the purpose of including thresholds. 
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To implement the super particle thresholds in the minimal low energy super gravity 
model, the renormalization group /^-function must be calculated in a form that allows every 
particle to be counted in the simplest possible way. This, of course, implies that we will 
have to make allowances for effective theories with half odd integer doublets, as discussed 
in Section VIII. We will, therefore, implement particle thresholds as steps in the particle 
content of the model. In the following, the one loop /3-functions of the gauge couplings and 
Yukawa couplings are considered. 



1. Gauge Couplings 

In the general case (but for a single, simple gauge group G), at one loop 

dg 1 



dt (47r) 

where 



, (Bl) 



h^"^ = \UF) + \US) - ^-^C,{G) , (B2) 

and where F, S", and G stand for the fermion, scalar, and adjoint representations, respec- 
tively. Using the fact that gauginos are in the adjoint representation gives in the S'[/(3) 
case 

= ^r2(F3) + ^r2(53) + ^C2(G3)% - yC72(G3) , (B3) 

where F3, 6*3 refer to the fermion and scalar representations, and = SU{3) and where 

r 1 , A > M, ; 

Mg is the mass of the gluino. In the SU (3) case, one has 

UR,) = 2(i)iVQ + (i)iV^+ (i)iV^ , (B5) 

where Np equals the number of families of particle p. This result is valid for both fermion 
{R = F) and scalar {R = S) representations. Equations (p3| ) and (p5|) lead to 

= liNu + N,) + ^Nq + ^N^ + liV^ + 2% - 11 , (B6) 

where we have assumed that Nq = {N^ + Nd)/2. Also the fact that left and right handed 
quarks of a given flavor have the same mass allowed us to use Nu = and = N^. Note 
that Eq. ( |B^ ) reduces to the right standard model result when Nq = = = 9g = 0, 
and to the right supersymmetric result (A > Msusy) when Nq = iV^ = At = 3 and 6g = 1. 
Similar formulas are calculated for gi and g2 
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^ 5^2 12 '' 4 ^ 4 30 ^15 "15 ^ 

+ 1^^^ + + J^^*" + ^^^^ + ]^^^*" + ^^'^^ ' ^^^^^ 

+^(iv^„ + iv^J + ^(iv*„ + iv*J + . (B7b) 

The two loop contributions may also be calculated in this manner. 



2. Yukawa Couplings 

We find it useful to define new doublet fields 

= sinP^u + cosP^d , (B8a) 
$iy = - cos - sin /3l>d , (B8b) 

where tan/3 = v^/vd and $ = iT2^*. The VEVs of these two new fields are < $° >= v and 
< $5/ 0- can now rewrite the Yukawa interaction Lagrangean of the MSSM in terms 
of these fields. The /3-function for the Yukawa couplings can be computed from formulas in 
Ref. [0. At one loop, the renormalization group equation for a general Yukawa coupling is 



given by 



dt {An) 
After replacing 



^ Y/?«. (B9) 



= sin/3l>fe - cos/?$H , (BlOa) 
= cos p<^h - sin (3^H , (BlOb) 



in the up sector Yukawa interaction Lagragian, Eq. ( |B9|) yields 

^ = ^(=YJA.v., . (BID 
where we will use the notation s = sin (3 and c = cos f3. From Eq. ( |B11| ), we get 

Using sin/5 = Vu/v and Eq. ( [A^ ) to one loop, we can write 

dins 2/^" 
—1— = c ( ) 

at Vu Vd 

^ c2(7$„-7$J. (B13) 



(47r 
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Putting Eqs. ( pil|) and (|B13| ) together yields 
dY,. 1 



dt (47r)2 

Similarly for the down and lepton sectors, we get 



[Y,AsYj-c^(7c.„-7*jYj . 



dY, 



dt (47r)2 
dY^ 1 



[Y,/3(eY,) + s'(7*„-7*jYj , 
[Y,/3(cYj + s2(7*„-7*jYj . 



(B14) 

(B15a) 
(B15b) 



dt (47r)2 

Now we list the forms of the /^-functions appearing in Eqs. (P14|) , (P315a|) , and (|B15b|) 



YeA 



\{,'e,, + c2^^)Y„YtY„ + \{,^e-, + c^^^)( 2{;^ y^yDy, 



k=l 



N- 



1 



+ YJJ: Yt YJ ) + -(c^^, + s^^^ - Ac\en - en))Y^Y\Y, 

k=l 



+ + s'^^h)YJE Y^Y,} + YJ (s^^^. + ^Oh) Tr{3YtYj 



fc=i 



+ o\e, - On) Tr{3YtY, + YIyJ ] - { ^^??[ ^ + 



'36 9 



4^ 



7V^ 



(B16a) 



hc^, + s2^^^)Y,YtY, + i(c2^^ + s^0^)( 2{;^ Y,Yt}Y, 
^ ^ fc=i 



At 



+ Y,{Y^ Y\Y,} ) + -(s^^, + c'Oh - As^On - ^^,))Y,Yt Y, 
k=i ^ 

1 

+ + c'^h)Y,{E YlYj + YJ (s2(^^, - ^^^) Tr{3Yt YJ 

^ fc=i 

+ (c^^, + s'Oh) Tr{3YtY, + Y^YJ ] " { ^ + ^^^^ 



(B16b) 



^ hc'e, + s^^^H)Y,YtY, + i(c^^^ + s'e^){ 2{Y^ Y,Yt}Y, 
^ ^ fc=i 



+ YJ5: YtYj ) + YJ {s\er. - en) Tr{3YtYj 



fc=i 
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+ + s'9h) Tr{3YtY, + YtYj ] - { -gl[ - + 



4 



} 



(B16c) 



where the various 6s equal zero below the mass threshold of the respective particle and equal 
one above it. Note that Eqs. ( [B16aj ), ( |B16b| ), and ( |B16c| ) reduce to Eqs. (|X9|), ( [A10| ), and 
( |A11 ), respectively, in the supersymmetric limit in which all ^s are equal to one. When 
6h = ^ and all other ^s in the above equations are equal to zero, we recover the standard 
model result after identifying the standard model Yukawas as Y^ 
Y": = cY.. 



sY., Y^ 



cY^, and 
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